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Abstract. The phonon spectral function of the one-dimensional Holstein model 
is obtained within weak and strong-coupling approximations based on analytical 
^H ' self-energy calculations. The characteristic excitations found in the limit of 

jy^ ' small charge-carrier density are related to the known (electronic) spectral 

properties of Holstein polarons such as the polaron band dispersion. Particular 
jrt , emphasis is laid on the different physics occurring in the adiabatic and anti- 

adiabatic regimes, respectively. Comparison is made with a cluster approach 
exploiting exact numerical results on small systems to yield an approximation 
for the thermodynamic limit. This method, similar to cluster perturbation 
theory, confirms the analytical findings, and yields accurate results also in the 
intermediate-coupling regime. 



PACS numbers: 71.38.-k, 63.20.Kr, 63.20.Dj, 71.38.Fp, 71.38.Ht 



Intermediate or strong electron-phonon (EP) interaction gives rise to the existence 



(N 

> _ 

^^ ' 1. Introduction 

o 

^s . of polaronic carriers in a number of interesting materials (see, e.g., [1-3]). As a 

\^ ' consequence, models for such quasiparticles, e.g., the Holstein model [4] considered 

^^ . here, have been investigated intensively in the past decades in order to understand 

the process of polaron formation. Whereas valuable insight into the latter may be 
^ , gained by considering a single charge carrier (i.e., the Holstein polaron problem), real 

y ' materials are usually characterized by finite carrier densities, motivating studies of 

many-polaron models [5-8]. 
^ ' Over the last decade, a large number of theoretical studies, the most reliable 

Q , of which based on unbiased numerical methods, have led to a fairly complete 

O ' understanding of the Holstein polaron concerning both ground-state and spectral 

properties. Calculations of the latter, e.g., the one-electron Green function, are 
particularly rewarding as they provide detailed insight into the non-linear process 
, . , of an electron becoming self-trapped in the surrounding lattice distortion. For reviews 

C^ ' of the Holstein polaron problem see [9, 10]. 

In this paper, we contribute to a completion of the knowledge about the 
single-polaron problem by investigating the phonon spectral function, an important 
observable rarely considered in previous work. On the contrary, it has been studied 
quite intensively for the spinless Holstein model and the Holstein-Hubbard model, 
both at half filling, for which phonon softening at the zone boundary has been found 
to occur at the Peierls transition in one dimension [11-16]. Note that the assumption 
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of a local self-energy frequently used in combination with dynamical mean field theory 
(DMFT) — as appropriate for infinite dimensions — leads to an unrealistic wavevector- 
independent softening of all phonon excitations [17,18]. The renormalization of 
the phonon modes in small clusters with one electron, especially their softening at 
the critical EP coupling, has been investigated numerically in [19], and results for 
the coherent phonon spectrum of the infinite system have been reported in [20]. 
Furthermore, analytical approximations for the phonon self-energy and the frequencies 
of the vibration modes in the coupled EP system have been obtained in [5, 21, 22]. 

Here we present analytical calculations valid at weak and strong EP coupling, 
respectively, as well as a cluster approach, similar to cluster perturbation theory [23- 
25], which yields accurate results in all relevant parameter regimes. Both approaches 
are capable of describing the momentum dependence of the phonon spectral functions, 
and we find a very good agreement between analytical and numerical results. 

Despite the formal restriction to the one-electron case, the numerical results will 
correspond to finite electron densities due to the finite underlying clusters. The 
analytical calculations presented here are based on the electron (polaron) spectral 
functions previously deduced in [26] for the weak-coupling (WC) and strong-coupling 
(SC) cases and for charge carrier concentrations n < 0.5. These spectral functions 
depend on n and their calculation requires a self-consistent determination of the 
chemical potential /i. To avoid this difficult task, we shall consider the analytical 
formulas in the mathematical limit of fj, approaching the bottom of the electron 
(polaron) band. The resulting limiting shape of the spectral functions will provide us 
with a picture of the phonon spectrum for small carrier concentrations. In particular, 
as demonstrated below, the positions of the low-energy excitations will be found to 
be in a very good agreement with numerical results. Moreover, the general analytical 
treatment will enable us to discuss to some extent the features which occur for non- 
negligible carrier concentrations. 

The paper is organized as follows. In section 2, we introduce the Holstein model. 
Section 3 is devoted to the derivation of the analytical results and a discussion of our 
cluster approach, whereas section 4 contains the results and a discussion. Finally, we 
summarize in section 5. 

2. Model 

In view of the SC approximation presented below, it is convenient to write the 
Hamiltonian of the one-dimensional (ID) spinless Holstein model as 

H = -t^Yl "^ + ^0 E ^h E ^rA^i ' (1) 

j j jj' 

where 

C,,^gu;o{b]+b^), C^,-,)=t. (2) 

In equation (1), c] (bl) creates a spinless fermion (a phonon of energy ujq) at site i, 
and fii = c^c^ with n^ = 0, 1. The first term contains the chemical potential /x and 
determines the carrier density n, whereas the second term accounts for the elastic and 
kinetic energy of the lattice. Finally, the last term describes the local coupling between 
the lattice displacement Xi ~ b\ + b^ and the electron density "hi with the coupling 
parameter g (for j = j' , cf equation (2)), as well as electron hopping processes between 
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neighbouring lattice sites (j'j) with hopping amphtudc t (for j = j' ± 1). We take t 
as the unit of energy, and set the lattice constant to unity. 

From numerous previous investigations [10] of this model in the one-electron case 
considered here the following picture emerges. At WC, the ground state consists 
of a large polaron, corresponding to a self-trapped electron with a lattice distortion 
extending over many lattice sites. As the EP coupling is increased, a cross-over takes 
place to a small-polaron state, in which the lattice distortion is essentially localized at 
the site of the electron, leading to a substantial increase of the quasiparticle's effective 
mass in the intermediate-coupling (IC) and SC regime. Depending on the value of 
the adiabaticity ratio a — cuo/t, the small-polaron cross-over occurs at a critical 
value of the EP coupling determined by the more restrictive of the two conditions 
A = Ep/2t > 1 (relevant for a < 1) or g^ = Ep/ujQ > 1 (for a > 1). Here Ep = g'^ujQ 
is the polaron binding energy in the atomic limit defined by i = 0. 

3. Methods 

3.1. Analytical approach 

The aim of the present treatment is to deduce and to interpret the essential features of 
the numerically calculated phonon spectral functions (see section 4). In the following 
calculations, we shall deal with coupled equations of motion of the Matsubara Green 
functions for the lattice-oscillator coordinates on the one hand, and for the spinless 
charge carriers on the other hand. It will be shown that the general features of the 
spectral functions may be understood on the basis of the results obtained for the WC 
and SC cases, where an approximate treatment is well justified [7, 26]. 

3.1.1. Weak-coupling approximation The Matsubara Green function for phonons, 
defined as 

D{miTi;m2T2) = -(Tt-Xtoi (ti)x,„2 (^2)) , (3) 

fulfills the equation of motion 
1 / 92 



o~ ( ^^ ^ ^0 ) -D(miTi; TO2T2) = (5mim2<5(Ti - T2) (4) 

- gUJo{TrC^^{Ti)cl^^{Ti)x,n2iT2)) , 

assuming {x.^^) = in the WC regime. The mixed Green function on the rhs of 
equation (4) will be expressed by means of the generalized fcrmionic Green function 

[27, 28] 

G(niTi;n2r2;;7) = -^(T,c„^(ri)4(r2)5)^G(l;2;;7) (5) 



and 



r^- CXp - / dr ^ Unn' {T)Cnn' (t) 



5 - T, CXp - / dr V Unn' {T)Cnn' (t) , (6) 



where the classical variables Unn' (t) were introduced as a purely formal device. 
Consequently, the relation between the mixed Green function and the fermionic Green 
function reads 

-G{niTi;n[T[;U) 



{TTCn,iTl)cl,iT[)Cm^m^{T2)) = 



SUm2m2{T2) 



(7) 
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Here S/SU denotes the functional derivative. Defining the inverse Green function 
G^^l; 2; U) to G(l; 2; U) according to [27] we obtain 

(^rC„^(Tl)cJ„^(Ti)C™,™2(T2)} (8) 

rl3 



f dr" VG(TOiTi;m"T") 
Jo , 



rP 

-G-\m"T";m'T';U) 



SUm2m2{T2) 



G{m'T';miTi) 
u=o 



for the interaction term in the equation of motion (4). The resulting equation for the 
phonon Green function is quite general, as no approximations have been made up to 
now. 

In the sequel, the fermionic Green functions in equation (8) will be obtained using 
the fermionic spectral functions A{k,uj) which have been calculated to second order 
in the EP coupling constant in [26]. To the same order, the functional derivative of 
G~^ will be determined using the relation between the inverse Green function and 
the self-energy E, G-^{1;2;U) = G^\l:2) - E(1;2;C/), where the Green function 
of the zeroth order. Go, is independent of U. Accordingly, the derivative of G^^ in 
equation (8) will be expressed as the derivative of S, which gives to second order 
[26,29,30] 

"' " " ' ' "~ =^ {TrCm"m'iT")G„i2rn2iT2)}Sm"m'S{T" - t') .(9) 

(7=0 



-^{m"T";m'T';U) 



5Um2m2{T2) 

Therefore, equation (4) acquires the following form: 



2uJo \dTl 



i^o] D{miTi;m2T2) ^ 5,nirn2S{Ti -T2) - dT'^G(miri;TO'T') (10) 

m' 
X {TrGrn'm'{T')Cm2m2{T2))G{m'T';miTi) . 

However, 

(7VG,„'„'(t')G„2„2(t2)} = g'^Ujl{TrXjn>{T')Xjn^{T2)) = - g^ UJ^ D {m' t' ; m2T2) . (11) 

Substitution of equation (11) into equation (10) and subsequent multiplication on the 
rhs by -D^^ give 

n(miTi; TO2T2) = .g^cjoG(TOiTi; m2T2)G(m2T2; miTi) (12) 

with n = Dq — I?~^, the self-energy of the phonons in the Matsubara framework. 
Fourier transformation of equation (12) leads to 

U{q,iujn) ==.g^Wo— ^-^G(fc,iw^)G(fc + (7,i(w^+w„)), (13) 

k '' V 

where ujy = {2v + l)7r//3 and cij„ = 2n7r//3 are Matsubara frequencies for fermions and 
bosons, respectively. 

The substitution of the spectral representation of the fermionic Green functions 

G(/fc',ic..)= r dc'^i^ (14) 
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into equation (13), the summation over the frequencies o;^, the analytical continuation 
iw„ >-^ Lu = Lu + i5 and the subsequent limit 5^0 lead to 

1 /'OQ /'OO 

X A(k + q, e%f{uj') - /(6')]C(^ + ^' - e') (15) 

with 

/(^) = ^^' ax) = --i7T5{x). (16) 

eP^ + 1 X 

In the low-temperature limit the integrand of equation (15) may be non-zero only if 

w' < and e' > ^ f{uj') - f{e') = 1 or 

tj' > and e' < ^ f{uj') - f{e) = -1 . (17) 

Using the second-order result for the electron spectral function A{k\ijj') deduced in 
[26], the calculation of the phonon self-energy (15), the corresponding retarded phonon 
Green function D^{q,Lo) and the phonon spectral function B{q,uj) would be straight 
forward [31]. However, A{k',LLi') as derived in [26] depends on the charge-carrier 
concentration n and has to be determined self-consistently with a condition for the 
chemical potential /i. The situation is thus simplified if we restrict ourselves to the 
case of small carrier concentration, for which the dependence on n is expected to be 
unimportant. To choose the dominant contributions to the integral on the rhs of 
equation (15) for small n, the integration over lo' and e' will be divided according to 
the character of the electronic spectral functions A{k, oj') and A{k + q, e'). 

The coherent part of the spectrum Ac(k,uj), non-zero for \uj\ < luq, consists of 
quasiparticle peaks, 

A,{k,u;)^ZkS{uj-{Ek-ti}), (18) 

whereas outside this frequency interval the incoherent spectral function Aic{k,ui) is 
formed by peaks of finite width. In the sequel, the integrals obtained in this way will 
be examined with respect to the behaviour in the limit of small concentrations, i.e., 
for fi lying near Eq, the bottom of the band defined by equation (18). 

The real and imaginary parts of n(g, oj) are obtained if the real and imaginary 
parts of the ^-function are substituted into equation (15). The behaviour of 
the resulting integrals for fi in the neighbourhood of Eq is then analyzed in the 
mathematical limit /i — > Eq . We find that the integration for Jicll(q,uj) according to 
equation (15) gives zero in the limit of vanishing carrier concentration. In contrast, the 
integration for Imn(g, lo) yields a non-zero result in this limit, with the coherent parts 
of the spectra, Ac{k, uj') and Ac{k + q, e'), being the only non- vanishing contributions. 
Taking into account only the latter, we shall express Imn((j, oj) as 

1 
T 



Imn(g, w) = - -g^ujl / dkzkZk+gd{uj - (Ek+q - Ek))9{n - Ek) 



LUq 



X / <i€'5{e' ^{Ek+q-^l)) 







-g ujo dkzkZk+q9{fi - Ek)-— ^^7- 



-g luq / dkzkZk+q9{n - Ek) 

'\k\<kr \^ko--^k„+q\ 

X I Ae' 5{e' - {Ek+q - p)) , (19) 
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with Ek < 11 < Ek+q < Wo- According to [26], the energies of the electronic band are 
given by the equation 

""'-^'^ ^^y_,, ^i-ii/wY E,-u,-^ ' ^^^^ 

where £,k = — VFcosfc [W — 2i), and its derivative 

£;^ = VFzfe sin fc , (21) 

with the spectral weight [26] 



Zk 



r?, i2 r^ AC 

.9 ^^ / Q? 



ttT^ ]-w ^l - i^/Wf {Ek -too- 0^ 



(22) 



The wavevector kg in equation (19) lies in the neighbourhood of fc = defined by 
9{n — Ek) ^ I and satisfies 

u; - Ek+g + Ek = 0. (23) 

Consequently, the value of ko appears to be a function of lu at fixed wavevector q, 
fulfilling Eko < M- 

To obtain a qualitative picture of the phonon spectral function, we explicitly take 
the limit /i -^ Eq . Thereby, lim ^^+ Kcll{q,u!) ~ and the rhs of equation (19) is 
non-zero only on the curve lu = LUq = Eg — Eq provided that < uj < ajQ. The values 
of liaIl{q,Lu) are given as follows: 

ImU{q,uj) = - 152^2 _io_^(^^ _ (^^ _ ii;o))A(w - {Eg - E,)) , (24) 

where the function A(xi — X2) is the generalization of the Kronecker symbol being equal 
to unity for xi = X2 and zero otherwise. Consequently, in this limit, the narrow peaks 
of finite frequency width following for Imn((7,a;) from equation (19), are replaced 
by discrete lines given by equation (24). The divergence of the imaginary part of 
the phonon self-energy at g = and g = tt is connected with the ID electron band 
dispersion. Our method of calculation only takes into account contributions to the 
electronic spectral function up to second order, which is insufficient if divergences 
occur. 

On physical grounds, one expects no renormalization of the phonon excitations 
in the zero-density limit, i.e., for one electron in an infinite lattice. The non-zero 
imaginary part of the phonon self-energy even for n — > [equation (24)] does not 
contradict this expectation since the integrated weight of the corresponding features 
is zero. As discussed in section 4, these non-zero contributions to the limit of the 
phonon spectrum may be related to results for small but finite carrier densities. 

So far, we have restricted ourselves to a; > 0, but the case w < may be treated 
quite analogously. The only non-zero contribution to lTaI{{q,Lo) in the limit ji -^ E^ 
is obtained for Ek+g < fx < Ek < ujq 'm the frequency range —ujq < o; < on the curve 
ujg = —{Eg — Eq), with the result 

ImU{q,oj) = -ImU{q,\uj\). (25) 

The retarded Green function D^{q,ui + iS) as the analytical continuation of 

D{q, iujn) = J. — TK 2 o — TTr~- — V (^6) 

(itj„)2 - t^if - 2tJon((7, iw„) 
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in the upper complex half-plane determines the phonon spectral function 

B{q,uj) = --lrRD^{q,Lj + iO+). (27) 

TT 

Using equation (26) and the preceding analysis of Il(q,u! + iO+) in the limit /i -^ Eq , 
we may conclude that for w > 

B(n , A - 1 (2^o)^Imn(g,^) 

if Lo ^ Eq ~ Eq and < uj < loq, where Imll is given by equation (24). Otherwise 
lmU{q,Lj > 0) = and 

B{q,uj)^S{uj-uJa). (29) 

The result for a; < 0, analogous to equations (28), (29) fulfills the relation 

B{q,Lj) = -B{q,-L,), (30) 

in agreement with the general requirement on the imaginary parts of retarded Green 
functions of real dynamical variables [32] . 

3.1.2. Strong- coupling approximation The equation of motion (4) does not hold 
exactly in the SC regime, as the coordinate Xm of the local oscillator at site m 
implies a shift due to the local lattice deformation associated with on-site small-polaron 
formation. In fact, {xm) = 2g{cl^c^), which follows from the Lang-Firsov canonical 
displacement transformation [33] . However, dealing with the limit of negligible charge- 
carrier concentration, (xm) = may again be assumed. On the other hand, the charge- 
carrier number operator in the electron picture is equal to the number operator in the 
small-polaron picture. Therefore, we interpret the Fermi operators c^j(ti), cI^_^{ti) 
in equation (4) as annihilation and creation operators of small polarons — the correct 
quasiparticles in the SC limit. Accordingly, the mixed term on the rhs of equation (4) 
will be expressed using the generalized small-polaron Green functions, defined again 
by equations (5), (6), where the operators Cnn(T) = gujQXn{T) as before, but the Ci^nn') 
correspond to the nearest-neighbour hopping term in the SC regime, i.e., 

q„„,) = iexp{-g(6t^ - 6„ - h^, + fe„,)} • (31) 

The formalism of generalized Green functions of small polarons was introduced by 
Schnakenberg [34] and applied to self-energy calculations in [29,30,35]. Apart from 
the C/„„/\ given by equation (31), in contrast to previous work we also include 
Cnn = g^o^n in our definition of the generalized Green function. 

The presence of the coefficients C„„ in the time-ordered exponential in 
equation (6) causes the polaronic operators not to commute with the exponent due 
to the oscillator shift proportional to cl^c^^. Accordingly, the zeroth-order generalized 
small-polaron Green function Go, corresponding to the atomic limit t = 0, is U- 
dependent because it fulfills the equation of motion 

rj- 2g'^ujoUnindTi) Go{niTi;n2T2;U) = (J„i„2(5(ti - T2) , (32) 



where r] — ^^ — g^^o- Using matrix notation [27], 



Go 1(1; 1'; (7) - f-A. - rj ~ 2g^uj„Un,nAri)] Sn,n[S{T, - t[) (33) 



represents the matrix inverse to Gq. 
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To obtain a qualitative picture of the phonon spectral function at SC we consider 
the limit g^ ^ 1. According to previous considerations [7,26], the polaron spectral 
function in this limit is dominated by the coherent part representing the polaron 
band of width 2We~^ , showing that small polarons are the correct quasiparticles 
and that multi-phonon processes in the self-energy are negligible. Consequently, the 
Green functions in equation (8) will be expressed by means of equation (14) using the 
coherent polaron spectral function 

A^{k,u) = 5{u-{^k+v)), (34) 

where ^fc = —W cosk {W — 2te^3 ) and 



^ G-\l-A'-U) 



0Um2m2V2) 



U=0 



0Um2m2 {'''2 



-Go-i(l;l';t/) 



(35) 

(7=0 



is assumed to be a good approximation. 

Using equation (35) and substituting equation (8) into equation (4), 

1 / a^ \ 

^ — ( ^"3 - ^0 ] D{miTi;m2T2) = <5„iim2(5(Ti - T2) (36) 

- 2g'^LUoG{miTi; m2T2)G{m2T2; miTi) 

is obtained. Fourier transformation of equation (36), use of the spectral representation 
of the polaron Green function based on equations (14), (34), and summation over the 
fermionic Matsubara frequencies result in 

2ujo 2 2(jJo 

(iuJn)'^ - UJq ' (it^n)^ - L^'q 



D{q, iw„) = - — ^2 2 ^ 2.9^^0- — -3 ^I{q, iu;„) (37) 



/(g, iiun) ^ 1 V /" d^' r de'A{k, u;')A{k + q, e'^ '^^'^'^ •^^'' 



N ^ J -co J -co ' ' tj' - e' + iw„ ■ 

The analytical continuation icun '-^Lu = Lu + i6 into the upper complex half-plane gives 
the retarded Green function D {q, ZU) and, according to the argument at the end of the 
preceding section, only a; > is to be considered. After the analytical continuation, 
the integral I{q,Lu + iO"'') on the rhs of equation (37) is analogous to the integral in 
equation (15), but the coherent spectrum (34) is limited to a frequency interval well 
below the value lo — luq- The equations determining Im I{q, uj) are simplified compared 
to equations (19) -(23) since, according to the small-polaron spectral function (34) 
used, the spectral weight Zk — 1 and the band energy Ek = S,k ^ .9^'^o- The limiting 
procedure described in section 3.1.1 gives again Re/(g, w) = and 

Im/(,,c.)^- ^^--;^^-f°» . (38) 

214^1 singi 

Consequently, the mathematical limit of Im D^{q,u!) for n -^ ^(0) — g'^ujo yields the 
limiting shape of the SC phonon spectral function as 

^/ ^ c-/ ^ 2g2tj2 A( LO - W (1 - COS q)) 

^0 ^ 7rT4^|sinq| 

The second term on the rhs of equation (39) reflects the energies of the small-polaron 
band which in the SC case lies entirely below the phonon resonance frequency luq- The 
divergences occurring at g = 0, tt are again related to the dispersion of the ID band 
and result from the failure of the approximation used at these wavevectors. 
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3.2. Numerical cluster approach 

The numerical approach used here is similar to cluster perturbation theory [23-25] for 
the one-electron Green function. For the case of the phonon Green function, it has 
first been proposed in [10] and applied to the half-filled, spinless Holstein model in 
[14]. A previous numerical cluster study of the renormalization of phonon excitations 
can be found in [19]. 

The phonon spectral function B{q,u}) is defined by means of the retarded 
phonon Green function which determines the response of the lattice to the external 
perturbation linearly coupled to the phonon variables [32]. The values of B{q, ui) may 
be shown to be proportional to the transition probabilities per unit time (at T ^ 
averaged with respect to the canonical distribution) for the transitions induced by the 
perturbation having frequency lu. 

In our case 

B{q,uj) = --lmD^{q,Lu) (40) 

TT 

will be calculated at zero temperature, so that the equality 

D^{q,Lj)^ \im{i;„\xq^-^ F7.^-,l^o) (41) 

r;^0+ U) + irj — II 

holds for uj > 0. Here \ipo) denotes the ground state of the infinite system, and the 
phonon coordinates are given by Xq — -^ ^ XjC^"'^. For the Holstein model (1) we 

have B{±q,uj) = B{q,Lj). 

The spectral function (40) fulfils a sum rule of the form J„ dujB(q,Lu) ~ 
(ipolXqX-qlipo), which relates the integrated spectral weight to the lattice elongation 
in the ground state [31]. The numerical techniques employed in the next sections 
guarantee that this sum rule holds up to machine precision. Note that the analytical 
results of section 3.1 fulfil the usual sum rule for the phonon spectral function [31] for 
T^O. 

To proceed, as a first step, we divide the infinite lattice into identical clusters of iVc 
sites each, and calculate the cluster Green function D^, [uj) of the Hamiltonian (1) 
with one electron and /i = for all non-equivalent pairs of sites i, j = 1, . . . , TVc- For 
this purpose, we employ the kernel polynomial method (KPM). Details about the 
computation of the Green function by the KPM and its advantages over the widely 
used Lanczos method can be found in [36]. The phonon Hilbert space is truncated 
[36] such that the resulting error of the spectra is negligible (< 10"'*), and we have 
used 1024 moments for the spectra shown below. 

In cluster perturbation theory, an approximation for G {k,io) of the infinite 
system is obtained by taking into account the first-order inter-cluster hopping 
processes, leading to a simple Dyson equation [24]. However, in the case of the phonon 
Green function, it turns out that the first-order term vanishes, since the electron 
number per cluster is conserved [14]. As a result, the cluster approach used here 
reduces to a Fourier transformation of the cluster Green function, 

i?^(q,^) = i-^Af^)(c.)e-^^(-^). (42) 

Nevertheless, it represents a systematic approximation to the exact Green function, 
as results improve with increasing cluster size Nc. Moreover, the method becomes 
exact both for a non- interacting system (g = 0), and in the atomic- or SC limit 
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t — 0. Wc shall see below that, provided N^. is large enough to capture the physically 
relevant non-local correlations, the method yields accurate results for all interesting 
parameters. Note that the defects mentioned in [14] originating from the neglect of 
true long-range order at half filling (above the Peierls transition) are absent in the 
low-density limit considered here. 

4. Results 

In contrast to previous considerations [5, 19-22] focussing on the renormalization of 
the vibration modes and their softening in the IC regime, the aim of this work is the 
calculation of the g-dependent phonon spectral functions in the entire frequency range 
and all relevant EP-coupling regimes. Particular attention is paid to the connection of 
the low-energy features with the electron (polaron) spectra studied previously [25, 26]. 
In the analytical calculations of section 3.1, the primary goal was the determination of 
the imaginary part of the phonon self-energy together with an analysis of the coherent 
and incoherent parts of the electron (polaron) spectral functions. 

For a clearer representation, in the figures, we shall only show the non-trivial 
lower excitations in the analytical results, i.e., B{q, uj < luq) according to equation (28) 
for WC, and only the second term on the rhs of equation (39) for SC, respectively. 
Furthermore, the analytical data will be rescaled for better visibility (see captions of 
figures 1 and 2). 

Weak coupling: If g = 0, only the resonance excitation of phonons having the 
eigenfrequency Uq = luq (the Holstein model (1) neglects any dispersion of the phonon 
branch) takes place, and the phonon spectral function is represented by the delta 
function (29). A non-zero EP coupling connects the lattice variables to the charge- 
carrier ones, giving rise to the low-frequency (off-resonance) part of the spectral 
function, which reflects the transitions to the excited polaronic states (of large or small 
polarons for WC or SC, respectively). According to the WC analytical calculations 
of section 3.1.1, this part of B{q,Lu) is given by equations (24), (28), and reflects the 
coherent part of the electron spectral function lying in the frequency range u> < uq. 
All this is confirmed in figure 1(a) for a = 0.4, i.e. in the adiabatic regime, which also 
shows the polaron band dispersion in the thermodynamic limit from variational exact 
diagonalization [37]. 

Contrary to the adiabatic case shown in figure 1(a), in the WC anti-adiabatic 
case (a = 4) reported in figure 2(a), the lower excitation in B(q,Lu) remains separated 
from the phonon line w = lvq, a-nd corresponds to the entire band of renormalized 
electron energies given — within our analytical approach — in Sec. 3.1.1 as the solution 
of equations (20) -(22). 

Both in the adiabatic (figure 1(a)) and the non-adiabatic (figure 2(a)) WC cases, 
we find a very good agreement of the WC approximation and the results from the 
cluster approach and exact diagonalization, respectively, with only minor deviations 
at large q for a = 4. These deviations, also affecting the g-dependence of the peak 
height in the analytical results, are a result of the shortcomings of the method for q 
lying near or tt (see section 3.1). 

An important point is that within the WC approximation, B{q,Lu) is strongly 
suppressed for q — and uj < uiq due to the divergence of equation (24). Hence, the 
peak in the numerical results is not reproduced. In contrast, the SC approximation 
corresponds to undamped quasiparticles (polarons) with a strong signal at g = (cf 
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equation (39)). Both these anomahes, connected with the dispersion in the ID electron 
(polaron) band, are a consequence of the approximations used and have no physical 
relevance. 

Intermediate coupling: The characteristic structure of B{q,Lu) consisting of the 
phonon line and the low-energy part continues to hold even at stronger EP coupling. 
Interestingly, for a — 0.4 at IC (figure 1(b)), we observe level repulsion between 
the weakly renormalizcd electron band and the bare phonon excitation at some wave 
number qy — determined by the point in q-space where the curves u = Eg — Eq and 
u = ojQ would intersect — as in the case of the electronic spectrum [25]. For A = 1 
(figure 1(c)), the critical coupling for small-polaron formation, the low-energy feature 
has already separated from the phonon line, the latter being overlaid by an excited 
"mirror band" lying an energy loo above the polaron band. 

The small-polaron cross-over for a = 4 is determined by the ratio g^ , and occurs 
at 5^ = 1. The phonon spectrum at this critical coupling is shown in figure 2(c). 
We detect a clear signature of the small-polaron band with renormalizcd half-width 
of about 0.70t, in good agreement with the SC result W « 0.74t, but an order of 
magnitude larger than in the adiabatic case of figure 1(c). In the latter, the SC 
approach predicts W w O.Oli, which is significantly smaller than the numerical result 
of about 0.09t. The fact that the analytical SC results are more accurate in the 
non-adiabatic than in the adiabatic IC regime has been pointed out before in [26]. 

Similar to figure 1(c), figure 2(c) features a mirror image of the lowest polaron 
band — shifted by loq — with extremely small spectral weight, which is barely visible for 
the number of Chebyshev moments — determining the energy resolution — used here. 

In the non-adiabatic case, small polarons exist even at intermediate EP coupling 
(see [7] and references therein). Therefore, as underlined by the above comparison 
of the bandwidths, the SC approximation agrees well with the exact results even for 
g < gc (figure 2(b)) and more than ever for g — Qc (figure 2(c)). 

By contrast, for a ~ 0.4, a large-polaron state exists for A « Ac, and both the 
WC (not shown) and the SC approximation fail to reproduce the characteristic IC 
features of figures 1 (b) and 1 (c) . 

Strong coupling: In the SC limit, the low-energy part is separated from the phonon 
line for both the adiabatic and the anti-adiabatic case (figures 1(d) and 2(d)), and 
the full spectral function has the form given by equation (39). Moreover, the effect of 
the polaron band-narrowing is more pronounced for the adiabatic case, as the small- 
polaron band at fixed A has the half-width W = We~'^^/'^° . In fact, there is no 
visible dispersion in the lower or upper band in figure 1(d). As expected, for these 
parameters, the SC approximation fits well the numerical results. 

Another general feature related to the WQ-dependence of the spectral function 
becomes apparent if we compare the heights of peaks in figures 1 and 2 (the ordinate 
scale in figure 1 is about a factor of five larger than that in figure 2). This dependence 
reflects the fact that the transition probability at low temperatures is proportional 
to the density of polaron states in the coherent band. The latter quantity increases 
with decreasing bandwidth, and the pronounced difference between adiabatic and anti- 
adiabatic spectral functions in the SC limit is evident from equation (39), where the 
height of the peaks is dominated by the factor 1/W ex c^^/'^". 

We conclude that the WC approximation based on the second-order electron 
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Figure 1. (colour online) Phonon spectral function B(q,w) as obtained by the 
cluster approach, for a = 0.4 and different EP couplings A and cluster sizes A^c- 

WC (a) and SC (b) - (d) analytical results ( , see text) have been multiplied 

by a factor of 10. Also shown are the bare phonon frequency ujq ( ), and the 

polaron band dispersion Eq — Eg ( — ■ — ) as calculated with the method of [37] . 
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Figure 2. (colour online) As in figure 1, but for a = 4 and for different couplings 

g^. WC (a) and SC (b) - (d) { , see text) analytical results are shown on a 

logarithmic scale. 
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spectral function, and the Hartree approximation for the SC small-polaron limit 
are able to grasp the main qualitative features of the phonon spectral function 
across the range of model parameters. However, as discussed in section 1, the 
numerical calculations do not represent exactly the limit of negligible charge-carrier 
concentration, owing to the restricted cluster volume. A direct comparison with 
results obtained by carrying out the limiting procedure for the analytical formulas 
in section 3.1 would be possible only for iV — > cxd. 

On the other hand, no restrictions concerning the charge carrier concentration 
were imposed in deducing equation (15) in section 3.1.1 and equation (37) in 
section 3.1.2. Consequently, the analysis of the latter equations outlined in section 3.1 
permits us to discuss the additional features of spectra for non-negligible carrier 
concentrations revealed by the numerical results. First, the discussion in section 3.1.1 
strongly suggests that the low-energy peaks of finite width correspond to the solutions 
of equations (19) - (23) at fixed wavevector g if /U lies above the bottom of the electron 
(polaron) band. Second, at finite concentrations, the contributions of the incoherent 
part Aic of the electron (polaron) spectral function to Imn(g, w) are not negligible and, 
in this way, additional maxima above u = loq occurring in the numerical results may 
be understood as originating from phonon- assisted processes implied in Aic- Finally, 
the non-zero Ilell{q,uj) in general causes a shift of the bare phonon line away from 
Lo — wq but, according to numerical results, the latter is not very pronounced in the 
WC and SC cases. 

5. Summary 

We have presented results for the phonon spectral function of the Holstein polaron in 
all relevant parameter regimes obtained by a reliable and systematic cluster approach 
similar to cluster perturbation theory. The characteristic features of the spectra have 
been discussed and successfully related to analytical self-energy calculations valid at 
weak and strong coupling, respectively. As far as a direct comparison is possible, our 
findings are in agreement with previous work on the phonon dynamics. 

In particular, we have pointed out the important differences between weak, 
intermediate and strong coupling, on the one hand, and between the adiabatic and the 
anti-adiabatic regime, on the other hand. As revealed by the analytical results, the 
phonon spectra of the Holstein polaron are dominated by the bare, unrenormalized 
phonon line and the renormalized polaron band dispersion. At intermediate coupling, 
additional features such as level repulsion and mirror images of the polaron band have 
been observed. Together with previous studies of the electron spectral function and 
the renormalization of phonon energies, this work provides a fairly complete picture of 
the spectral properties of the one-dimensional Holstein polaron, which has been in the 
focus of intensive investigations over several decades due to the wide-spread relevance 
of polaron physics. 
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